Negative thermal expansion in NbF_3 and NbOF_2:A comparative theoretical study

Thermal expansion control is always an obstructive factor and challenging in high precision engineering field. Here,the negative thermal expansion of Nb F_3 and Nb OF_2 was predicted by first-principles calculation with density functional theory and the quasi-harmonic approximation(QHA). We studied the total charge density, thermal vibration, and lattice dynamic to investigate the thermal expansion mechanism. We found that the presence of O induced the relatively strong covalent bond in Nb OF_2, thus weakening the transverse vibration of F and O in Nb OF_2, compared with the case of Nb F_3.In this study, we proposed a way to tailor negative thermal expansion of metal fluorides by introducing the oxygen atoms.The present work not only predicts two NTE compounds, but also provides an insight on thermal expansion control by designing chemical bond type.     

求解矩阵方程AXB+CXD=F参数迭代法的最优参数分析

本文针对求矩阵方程AXB+CXD=F唯一解的参数迭代法，分析当矩阵A，B，C，D均是Hermite正（负）定矩阵时，迭代矩阵的特征值表达式，给出了最优参数的确定方法，并提出了相应的加速算法.     

Shedding vortex simulation method based on viscous compensation technology research

Owing to the influence of the viscosity of the flow field, the strength of the shedding vortex decreases gradually in the process of backward propagation. Large-scale vortexes constantly break up, forming smaller vortexes. In engineering, when numerical simulation of vortex evolution process is carried out, a large grid is needed to be arranged in the area of outflow field far from the boundary layer in order to ensure the calculation efficiency. As a result, small scale vortexes at the far end of the flow field cannot be captured by the sparse grid in this region, resulting in the dissipation or even disappearance of vortexes. In this paper, the effect of grid scale is quantified and compared with the viscous effect through theoretical derivation. The theoretical relationship between the mesh viscosity and the original viscosity of the flow field is established, and the viscosity term in the turbulence model is modified. This method proves to be able to effectively improve the intensity of small-scale shedding vortexes at the far end of the flow field under the condition of sparse grid. The error between the simulation results and the results obtained by using fine mesh is greatly reduced, the calculation time is shortened, and the high-precision and efficient simulation of the flow field is realized.     

水化硅酸钙纳米压痕分子动力学模拟方法优化

针对水化硅酸钙纳米压痕模型忽视了压头与基底之间相互作用的问题，由尺寸差异引起的金刚石压头难以计算的问题，以及Wittmann模型无法得到实际接触面积的问题，提出了新的模型与计算方法．结合分子动力学方法，采用金刚石压头-Wittmann模型基底的组合方式构建无定形态水化硅酸钙纳米压痕试验模型．在建模阶段，考虑到压头模型与基底模型粒子间尺寸差异，提出了等比例替换模型，通过公式推导并就不同尺寸模拟结果验证了等比例替换模型的可行性．在计算阶段，提出了局部前处理的弛豫方法进行模拟．确定最大荷载位置处的接触面积为546 nm2，进而求出水化硅酸钙模型硬度H为0.84 GPa、折合模量Er为30.52 GPa．并通过纳米压痕试验，验证了模拟结果的准确性，证明了模型的科学性，对今后水化硅酸钙(C-S-H)纳米层面的模拟具有重要借鉴意义．     

移动加热器法生长碲锌镉晶体的组分输运与界面形貌研究

移动加热器法(THM)生长碲锌镉晶体时,界面稳定性对晶体生长的质量有很大影响。本文基于多物理场有限元仿真软件Comsol建立了THM生长碲锌镉晶体的数值模拟模型,讨论了Te边界层与组分过冷区之间的关系,对不同生长阶段的物理场、Te边界层与组分过冷区进行仿真研究,最后讨论了微重力对物理场分布的影响,并对比了微重力与正常重力下的生长界面形貌。模拟结果表明,Te边界层与组分过冷区的分布趋势是一致的,在不同生长阶段,流场中次生涡旋的位置会发生移动,从而导致生长界面的形貌随着生长的进行发生变化,同时微重力条件下形成的生长界面形貌最有利于单晶生长。因此,在晶体生长的中前期,对次生涡旋位置的控制和对组分过冷的削弱,是THM生长高质量晶体的有效方案。     

Ab initio and Monte Carlo studies of phase transitions and magnetic properties of YCrO3: Heisenberg model

By using the density functional theory (DFT) and Monte Carlo simulations (MCS) with the Heisenberg model, we have studied magnetic properties of the bulk perovskite YCrO₃. The exchange couplings of the Heisenberg model and the magnetic anisotropy are investigated. The 110 direction in the crystalline structure of the compound has shown the minimum energy, it is the easy magnetic direction. Using Monte Carlo simulations, the magnetizations behavior, the effects of system parameters and the critical exponents of the compound YCrO₃ are implemented. It is shown that the bulk perovskite YCrO₃ belongs to the 3D Heisenberg universality class.     

基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.     

A Theoretical study on the charge and discharge states of Na-ion battery cathode material,Na1+xFePO4F

Na₂FePO₄F is a promising cathode material for a Na-ion battery because of its high electronic capacity and good cycle performance. In this work, first principle calculations combined with cluster expansion and the Monte Carlo method have been applied to analyze the charge and discharge processes of Na₂FePO₄F by examining the voltage curve and the phase diagram. As a result of the density functional theory calculation and experimental verification with structural analysis, we found that the most stable structure of Na_1.5FePO₄F has the P2₁/b11 space group, which has not been reported to date. The estimated voltage curve has two clear plateaus caused by the two-phase structure composed of P2₁/b11 Na_1.5FePO₄F and Pbcn Na₂FePO₄F or Na₁FePO₄F and separated along the c-axis direction. The phase diagram shows the stability of the phase-separated structure. Considering that Na₂FePO₄F has diffusion paths in the a- and c-axis directions, Na₂FePO₄F has both innerphase and interphase diffusion paths. We suggest that the stable two-phase structure and the diffusion paths to both the innerphase and interphases are a key for the very clear plateau. We challenge to simulate a nonequilibrium state at high rate discharge with high temperature by introducing a coordinate-dependent chemical potential. The simulation shows agreement with the experimental discharge curve on the disappearance of the two plateaus. © 2018 Wiley Periodicals, Inc.     

Novel local smoothness indicators for improving the third‐order WENO scheme

The local smoothness indicators play an important role in the performance of a weighted essentially nonoscillatory (WENO) scheme. Due to having only 2 points available on each substencil, the local smoothness indicators calculated by conventional methods make the third‐order WENO scheme too dissipative. In this paper, we propose a different method to calculate the indicators by using all the 3 points on the global stencil of the third‐order WENO scheme. The numerical results demonstrate that the WENO scheme with the new indicators has less dissipation and better resolution than the conventional third‐order WENO scheme of Jiang and Shu for both smooth and discontinuous solutions.     

Domain decomposition methods in a geometrical multiscale domain using finite volume schemes

The heat equation is solved by using a finite volume discretization in a domain that consists of a two-dimensional central node and several one-dimensional outgoing branches. Several interface connection options to match the submodels set on the node and on the branches, with or without continuity, are looked at. For each of them, a monolithic scheme is defined, and existence and uniqueness of the solution is proved. New schemes are deduced, which are obtained through domain decomposition methods in the form of interface systems, with one or two unknowns per interface. A comparative systematic study is carried out from an algebraic and numerical point of view according to the interface conditions: Dirichlet, Neumann, or Robin. An efficient diagonal preconditioning is proposed.